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The superfluid-Mott-insulator phase transition of ultracold spin-1 bosons with ferromagnetic and antiferro- 
magnetic interactions in an optical lattice is theoretically investigated. Two counterpropagating linearly polar- 
ized laser beams with the angle 6 between the polarization vectors (lin-S-lin configuration), driving an Fg = 1 
to Fe = 1 internal atomic transition, create the optical lattice and at the same time couple atomic ground states 
with magnetic quantum numbers m = ±1. Due to the coupling the system can be described as a two-component 
one. At ^ = the system has a continuous isospin symmetry, which can be spontaneously broken, thereby fix- 
ing the number of particles in the atomic components. The phase diagram of the system and the spectrum of 
collective excitations, which are density waves and isospin waves, are worked out. In the case of ferromagnetic 
interactions, the superfluid-Mott-insulator phase transition is always second order, but in the case of antiferro- 
magnetic interactions for some values of system parameters it is first order and the superfluid and Mott phases 
can coexist. Varying the angle 9 one can control the populations of atomic components and continuously turn 
on and tune their asymmetry. 

PACS numbers: 03.75.Lm,03.75.Mn,71.35.Lk 



I. INTRODUCTION 

The investigation of quantum phase transitions (QPT's) 
in optical lattices has become one of the major subjects in 
the physics of ultracold atoms, especially after a convinc- 
ing experimental realization in the system of spin-polarized 
^^Rb 1 1|. One of the interesting directions in this rapidly de- 
veloping field is the study of QPT's in ultracold bosonic atoms 
with spin degrees of freedom. Systematic studies of QPT's 
in spin- 1 bosonic systems with antiferromagnetic interactions 
have been performed very recently. It has been shown that 
there are different Mott insulating phases such as spin singlet, 
nematic, and dimerized ones \2se ^ Ai ^ ^nd the mean- 
field analysis shows that the transition between the singlet 
and nematic phases is weakly first order Fragmented 
condensates with topological excitations have been studied in 
Ref. H. In the presence of a uniform magnetic field, different 
polar superfluids with spins aligned along or opposite to the 
field direction are possible |7]. QPT's of spin-2 bosons have 
been studied as well . 

In the present paper, we shall investigate QPT's of spin-1 
bosons with coupled ground states in an optical lattice and 
their collective excitations. We consider a generalization of 
the setup we have studied recently |9], where the population 
of only two degenerate ground states with magnetic quantum 
numbers m — ±.1 was considered. The presence of atoms in 
all the three ground states leaves two low-lying linear com- 
binations, which transform in the fundamental representation 
of an isospin SU(2) which determines the symmetry and the 
dynamics of the resulting system. Varying the laser polar- 
ization one can continuously turn on and tune the asymmetry 
between the components and control the populations of com- 
ponents in the superfluid and the Mott phases. As we shall 
show, a rather rich scenario of quantum phase transitions of 
first and/or second order becomes accessible to experiments 
in this way. In the present paper, we explore the phase dia- 
gram and the low-lying excitations of the various phases. In 
Section II we introduce the laser configuration and the Hamil- 



tonian which describes the system. The Bose-Hubbard model 
for this system is obtained in Section III, and its isospin sym- 
metry is discussed. Section IV gives a derivation of the phase 
diagram in the mean-field approximation and contains our re- 
sults on possible coexistence of Mott and superfluid phases 
and first- and second-order phase transitions. In Section V 
the collective excitations consisting of Bogoliubov mode and 
isospin waves with or without a gap are analyzed. The final 
section summarizes our conclusions. 



II. SYSTEM 

We consider a system of neutral polarizable spin-1 bosonic 
atoms of mass M, possessing three Zeeman-degenerate inter- 
nal ground and excited electronic states characterized by the 
magnetic quantum number to = 0, ±1 {Fg ~ = 1), in 
a (quasi-) one-dimensional optical lattice. The latter is as- 
sumed to be created by two counterpropagating linearly po- 
larized laser waves of equal amplitudes and frequencies with 
the wave number k^, and the angle 9 (with < 9 < tt/2) be- 
tween the polarization vectors (lin-(?-lin configuration) 
and with detuning A from the internal atomic transition. The 
running laser waves form left- and right- polarized standing 
waves with the Rabi frequencies 

n±{z) = nocos{kLZ±9/2) , (1) 

which couple internal ground and excited states by V and A 
transitions. In order to avoid decoherence due to spontaneous 
emission, | A | must be much larger than the spontaneous emis- 
sion rate. Because of the large detuning the excited state can 
be adiabatically eliminated. The resulting effective Hamilto- 
nian for the atoms in the three degenerate ground states has 
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the following form ifTllfll : 



iJ = 



2Af az 'A 



(2) 



where V'q(z) is the bosonic field annihilation operator for the 
atom in the hyperfine ground state \F — l,a). Summation 
over repeated spin indices a, is implied. ^ajS is the vector 
of F = 1 (3 X 3) matrices in the representation furnished by 
the three degenerate ground states. The matrix 



n+n^ 



nl 

n,n_ 



where 



fit = (1 + cos 6* COS 2kLz) 



(3) 



(4) 



determines the lattice potential with the period Tr/k^ and at 
the same time couples the atomic ground states with m = 
±1. The parameters gs.a describe the repulsive interaction 
of the condensate atoms and the spin-changing collisions. 
In a (quasi-) one-dimensional case they have the form lUsIl 
gs,a — ^fi'as,al Ma\, where a^^a are symmetric and anti- 
symmetric scattering lengths, and — \j2h/Mu)i_ is the 
size of the ground state for the harmonic potential with the 
frequency lox_ confining the atomic system in the transverse 
directions. In our previous work |9], we have studied the case 
when the ground state with m = is not occupied. In the 
present paper, we shall deal with the situations when this state 
is occupied as well. 



III. BOSE-HUBBARD MODEL 

The V and A laser-induced transitions lead to two sets of 
orthogonal Bloch eigenmodes which we denote by the indices 
and A, respectively. We assume that the atoms stay always 
in the lowest Bloch bands with the dispersion relations EQ{k) 
and Ef^{k). As explained below, in the present work we are 
interested in the case of A < and w 0. In this special case 
the lowest-band approximation is well justified, because the 
lowest Bloch bands are separated from the first excited ones 



by an energy of the order of y i?fl?ir2Q/| A|, where Ej^ = 

hujji is the recoil energy. Then the spinor-field operator * can 
be decomposed in the Wannier basis localized at the minima 
of the lattice potential labelled by i. 



(5) 



=0,A 



where Wcri(z) = Wct(z — z.^) are three-component Wannier 
spinors for the lowest energy bands. They are obtained by the 
solution of the eigenvalue problem for the Hamiltonian (|2ji in 



the case gs^a = and satisfy the orthonormality condition 

/ W^j(z) • Wcr'j{z) dz — 6ijSa-a'- They have the form 



Wo. = (0,Wo„Of , Wa, = (W^+„0,Ty_,)' 



(6) 



The indices i, j label the sites of the one-dimensional periodic 
lattice. The Ua-i is the Bose annihilation operator attached to 
the ith lattice site. The phases (p^i are not yet defined and their 
proper choice is discussed in a moment. The hyperfine spin 
variable at lattice site i is given by 



with 



^aa'^= j dz (z) • F, • W,,, (z) 



(7) 



(8) 



Substituting Eq. (|5jl into Eq. (|3 and taking into account 
only the hopping between the nearest lattice sites and the on- 
site atomic interactions, we obtain the two-component Bose- 
Hubbard Hamiltonian 



<2J> 



X ^^ncnincri - 1) + K^no, 



i 



(9) 



where is a chemical potential. The tunneling matrix ele- 
ments 



wL+i(^)- 



2Mdz^ A 
the atomic interaction parameters 



1 •w,,,(z)dz , 

(10) 



Ua 

Uo 
K 
P 



^ dz , 



(.gs +5a) ( 

gs J iWo^l^ dz , 

{9s+ga) J \Wo^\^ {\W+,\' 

^9a I {W^,fW+,W^,dz , 



(11) 



dz 



and the relative shift of the mean energies of the eigenmodes 

6 = I Wi iz) . i-^^, + ■ Wa,(z) dz 



Wj,,(z) 



2M dz^ 

' 2M dz^ 



= — / [EA{k)~Eo{k)] dk 
Jo 



Wo,j(z)dz 

(12) 
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can be simultaneously changed by varying the laser intensity 
(which is proportional to CIq) and/or the angle 6, but the vari- 
ations of Jo- and 6 are much faster. The parameter P can be 
either positive or negative depending on the sign of the anti- 
symmetric coupling ga- 

The phases (pa-i are determined from the requirement of 
minimal energy of the Bose-Hubbard Hamiltonian, which 
amounts to demanding that J^, cos((^c,i — 'Pa,i+i) and 
P cos{2ipQi — 2ipAi) are maximal. This requirement leading 
to 



ipoi = 



f (f^i if > 
\(f„i±Tr if < 

if < 

LpAi ±tt/2 if ga> 



(13) 
(14) 



has been taken into account in the derivation of the Hamilto- 
nian (|9j, which therefore only depends on the absolute values 
I Jct| and \P\ besides the other system parameters Ua, which 
are positve in the examples we consider and of about equal 
size, K, which can be larger or smaller than the Ua depending 
on the sign of the antisymmetric coupling ga, and S. As one 
can clearly see from Eq. Jq depends on f2§ cos 6'/| A| and it 
is always positive 1 14]. Innrinciple Ja can also become neg- 
ative for large enough 9 but at such values of 9 the shift 
S becomes comparable to the atomic interaction parameters 
in Hamiltonian (|9jl and it becomes energetically more favor- 
able for atoms to stay in the mode. In the present work, we 
are interested in the situations when both modes are occupied. 
Therefore we have to restrict ourselves to small values of 9 
and consider Ja as well as Jq as positive quantities. Thus in 
the following always the first of the alternatives in Eq. iVH 
applies. In general, Jq > J\. 

In the present paper, we investigate only the case of red 
detuning (A < 0), since in this case the minima of the ef- 
fective periodic potentials for the A and modes coincide 
which leads to higher values of K and \P\ compared to the 
case of blue detuning (A > 0). In addition if 6' = 0, 
W+i = = Woi/V2 and we have 

Jo = Ja = J , Uo = Ua = U , K = U + P , S = 0. (15) 

The and A modes are degenerate in this case, because 
Eo{k) = E\{k) for any value of the laser intensity. Indeed 
the Hamiltonian then has an enhanced U(l)xU(l) symme- 
try. This can be seen by defining the global generators of an 
SU(2) algebra (with L the number of lattice sites) 



To = 



1 ^ 

L 

iz {^Ai^o^ - a^aA*) , (16) 
1 ^ 



IaM, 



In terms of these generators the density of the hyperfine spin 
variable f — J^i ^i/L for 9 = 0, which always points in the 
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FIG. 1: Dependences of the parameters in the Bose-Hubbard Hamil- 
tonian on e for *''Rb, q = CiI/AlorA = -5. 



X direction defined by the direction of the linear polarization 
of the external laser, is given by f = 2Ti{^2)^x/ L for ga < 

(9a > 0). 

Together with the number density h = X^i^Li ('^oi'^oi + 



aAi)/L, the T generates an U(2) algebra. However, for 
ga > {ga < 0) the Hamiltonian (|9} only commutes with the 
generators n and Ti (T2) of the respective U(l)xU(l) sub- 
groups which they generate. 

Increasing the angle 6, one changes the form of disper- 
sions Ecr{k) and Wannier spinors Wa-i{z), and these changes 
are different for the A and modes. As a consequence, the 
values of the parameters in Hamiltonian (|9} will be modi- 
fied (see Ref. 1 14| for Jo 



'^i^a/Aj cos 61 , Ref. |9| for Ja(6'), 

and Fig.^for the others) and relations dlSI l become violated. 
Therefore, varying the angle 9 one can continuously tune the 
parameters in the Hamiltonian (|9} and remove the degeneracy 
by reducing the symmetry to the U(l) generated by n. 

Throughout the paper in our numerical estimations we have 
used known values of scattering lengths of ®^Rb (ga < 0) and 
^^Na (ga > 0) lEl and uj± = 2tt x 200 Hz. The laser wave- 
lengths Al were chosen to be 780 nm and 600 nm, respec- 
tively, which corresponds to the Dl line. 

Despite the fact that Hamiltonian (|9} has been explicitly 
derived here only for a one-dimensional lattice, it can be real- 
ized in higher-dimensional optical lattices as well. This can be 
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done with the aid of additional pairs of Hnearly polarized laser 
beams propagating along the x and y axes. In order to avoid 
the interference of laser beams propagating in the orthogonal 
directions their frequencies must be sufficiently different. In 
order not to change the number of components in the Bose- 
Hubbard Hamiltonian, the scheme of the internal atomic tran- 
sitions should remain the same. This can be achived if there is 
no TT component of the laser field, i.e., all the laser beams must 
be polarized in a plane perpendicular to the quantization axis, 
which is chosen to be the z axis. Therefore the polarization 
vectors of the laser beams propagating along the y and x axes 
are to be parallel to the x and y axes, respectively. As dis- 
cussed above and confirmed by numerical calculations in the 
next section, the angle 9 must be very small and then Ja ~ Jo ■ 
Thus it is not difficult to adjust the laser intensities and detun- 
ings of the additional laser beams such that the anisotropy of 
the tunneling matrix elements in the higher-dimensional ver- 
sion of the Hamiltonian (|9jl will be negUgible. 




, - 

ffrA 0.25 



FIG. 2: Ground-state energy of the Hamiltonian <18t for *^Rb (ga < 
0) in the case 9 ^ 0. q ^ Q,'^/4ujrA = -4.5, fi/U = 0.5 (a), 
H/U = 2.5 (b). 



superpositions of the number states 



N N 

E E 

no— nj^—Q 



(19) 



IV. MEAN-FIELD APPROXIMATION AND THE 
GROUND-STATE ENERGY 



We have done calculations with N 
discussed below. 



10 and the results are 



In order to investigate the superfluid-Mott-insulator phase 
transition, we employ the mean-field approximation 1 15, l3l 

aLa^Tj ~ (a^j + al.'j - , (17) 

where ^/j^ is the order parameter for Bose-Einstein condensa- 
tion in the component a. The tp^ carry phases, which are in 
equilibrium either the same for the both components or differ 
by TT, as can be easily shown by energy minimization. Thus 
we may take the i/^cr as real without unwarranted loss of gen- 
erality. In this approximation, the Bose-Hubbard Hamiltonian 
becomes local and every lattice site is described by the Hamil- 
tonian 

H'BH = -^dY,JA{^l + ^-)^'--^l] (18) 

CT 

+ ^ ^n„{n„ - 1) + KhofiA 

a 

-■Li (^aldlcLAaA + dj^aj^aoflo^ - Sho - fi'^fia- , 

a 

where d is the dimensionality of the lattice and we have omit- 
ted the site index i. Since in the mean-field approximation d 
plays a role of a numerical parameter, we have done the cal- 
culations presented below for d = 1. In the case of d = 2, 3, 
there will be only numerical corrections, but all the qualitative 
features will remain the same. 

The values of the order parameters minimizing the ground- 
state energy Eg{ipo,tpA) of the Hamiltonian (I18t determine 
the phase diagram of the system at zero temperature: i/^cr = 
{ipa- 7^ 0) corresponds to the non-condensed (superfluid) 
phase of the cr-component. 

The ground-state energy i?g(?/'o, ■0a) h^s been calculated 
by numerical diagonalization of the Hamiltonian ( I18> using 



A. Ferromagnetic interactions (ga < 0) 

Typical dependences of the ground-state energy of the 
Hamiltonian ( I18> on the order parameters ipo, ipA for ga < 
are shown in Figs. |2] and |3l If 6* 0, the ground-state energy 
Eg is independent of the sign of the real 4'a, 

Eg{i;A,^o)^Eg{\^Al\M) . (20) 

and it is a symmetric function of "00 and ipA- Eg {ipo, iPa) — 
Eg {tpA , "00 ) • More generally, for complex tpcr the isospin sym- 
metry at 6* = demands (still for ga < 0) that Eg only de- 
pends on the invariants |0oP + iV'Ap ™d (V'aV'o + '^'oV'a)- 
One has to distinguish the following two cases, (i) There is 
a single minimum at tpo = 4'A = [Mott phase. Fig. |2ja), 
corresponding to a state of unbroken U(l)xU(l) symmetry], 
(ii) There are four equal minima at = ?/;^ 7^ [symmetric 
two-component Bose-Einstein condensate (BEC), Fig.lJJb)]. 
In any of these minima the system has a nonvanishing x com- 
ponent of the expectation value of the hyperfine spin density 

^ {T,) ^ irAi'o + roi'A)/L. 




FIG. 3: Ground-state energy of the Hamiltonian <18> for *^Rb (ga < 
0) in the case ^ = 1°, g = ^I/AurA = -4.5, ii/Ua = 1.71. 
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If 7^ the scales of ipo and ip\ become different. Eg is not 
a symmetric function of ■(/'o and -i/'A anymore, but Eq. (I20> re- 
mains valid. The case (i) is possible again. Due to the different 
scalings of -00 and ipA instead of (ii) we have ipQ ip\ 
(two-component BEC). In addition, if 7^ 0, there are ad- 
ditional cases with two equal minima of Eg reached at (iii) 
ijjA = 0, V'o 7^ [mixture of noncondensed component A 
and BEC of component 0, Fig. 13 ■ As we see in Fig.|4] the 
mean occupation numbers of the A component turn out to be 
so small in these regions that there seems to be no hope to 
observe the noncondensed component experimentally. 

The phase diagram in the case 6* = is shown in Fig. |5] 
The lobes correspond to the Mott phases with different parti- 
cle numbers per lattice site, and outside the lobes we have a 
superfluid phase. The superfluid-insulator transition is always 
second order for < 0. 



B. Antiferromagnetic interactions (ga > 0) 

In the case of positive ga, i.e., positive P, the situation is 
quite different. If 6* = 0, the only condition in the ground 
state is V'oV'A — V'aV'o = 0, which restricts the expectation 
value of the isospin (T) to the (1, 3) plane. Apart from this 
restriction Eg is a rotationally symmetric function of the or- 
der parameters for any values of /i and the dependences de- 



li/U 
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FIG. 5; Phase diagram for * Rb (ga < 0) in the case 61 = 0. 



picted in Fig. |2jb) do not occur. Case (i) of the previous 
subsection remains possible, but instead of (ii) the minima 
of Eg are located on a circle V'o + V'a = V'^ on which a 
point is selected by spontaneous symmetry breaking [stable 
two-component Bose-Einstein condensed phase with sponta- 
neously broken isospin symmetry with (T) oriented arbitrar- 
ily in the (1, 3) plane, Fig.|6jc)]. In addition, there can be not 
only global minima, but also local ones either at ■0a = V'o = 
[metastable Mott phase, Fig.|6jb)] or on a circle [metastable 
superfluid phase, Fig.|6ja)], depending on the system parame- 
ters. This means that in some parameter regions the superfluid 
and Mott phases can coexist and the phase transition can be 
either first or second order in contrast to the case of ferromag- 




FIG. 6: Ground-state energy of the Hamiltonian <18> for ^^Na (ga > 
0) in the case 6 = 0. q = Q'^/AujrA = -6.1, ii/U = 1.57 (a), 
H/U = 1.6 (b), n/U = 2.2 (c). 
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FIG. 7: Ground-state energy of the Hamiltonian <18t for ^^Na (ga > 
0) in the case e = 0.1°. q = ng/4LjflA = -6.1, fi/UA = 1.6. 




FIG. 9: Phase diagram for Na (ga > 0) in the case ^ = 0. The re- 
gions of metastable superfluid phase coexisting with the stable Mott 
phase and that of matastable Mott phase coexisting with the stable 
superfluid phase are denoted by MSF and MM, respectively. 



C. Mott phase, arbitrary ga 



netic interactions, where it is always second order. 

If 6* 7^ 0, the symmetry between the two components is 
destroyed and the scalings for ■4>q and V^a become different. 
The minima of Eg can be reached at either tpo = ipA = Q 
or ipA — 0, ipo ^ Q as in the case ga < discussed above. 
In addition, there are dependences shown in Fig. with (or 
without) one local minimum at -0a = V'o = and two global 
minima at = 0, 7^ 0. In contrast to the case of negative 
Qa, the populations of the components can be comparable in 
all the cases (Fig.|8}. 

The phase diagram in the case 6* = is shown in Fig.|9] It 
consists of a series of lobes corresponding to the Mott phase 
as in the case of ferromagnetic interactions. However, the 
superfluid-insulator transition for ?i = 2 is now first order, 
while for n — 1, 3 it remains to be second order In the case 
n = 2, the stable (metastable) superfluid phase coexists with 
metastable (stable) Mott phase in a small region near the phase 
boundary. 
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FIG. 8: Extremal values of the order parameter t/jo and the mean 
occupation numbers of A and components for ^^Na {ga > 0). 
q = -6.5, e = 0.1°. 



In the Mott phase, ipA — 'fpo = and for every number 
of atoms per site n there are n + 1 eigenstates of the local 
Hamiltonian dlSl l. For n < 3 the expressions for eigenstates 
are very simple. If the number of particles per site n is less 
than two, the terms which contain |P| in Eq. ( I18> do not give 
any contribution and we have the usual Fock states, 

^0 = 0, |n = 0) = |0)a|0)o , 
Ei,o = -i^^ + 6) , |n = l)o = |0)a|1)o , 
Ei,i - |n-l>i - |1)a|0)o . (21) 

The lowest of these energies gives the ground-state energy 
per lattice site. The states with higher energies correspond to 
excited states whose positive excitation energies are the dif- 
ference of the respective energies E and that of the ground 
state. In the case 6 ^ 0, S vanishes. The states \n — 1)q and 
\n = l)i become degenerate and any linear combination of 
these states forms the SU(2)-symmetric manifold of degener- 
ate ground states in the case of one atom per lattice site. This 
case requires a positive chemical potential > 0, and, as we 
shall see by comparing with the energy levels for the case of 
two atoms per lattice site, /i < U — \P\. Thus, in the Mott 
phase with one atom per site, the mean numbers of atoms in A 
and modes tt-a and uq are undefined, but the total number of 
atoms n = riA + uq is fixed [see Fig.l^a)]. If 6 ^ 0, the de- 
generacy is removed and n\ — 0, uq — 1 [Figs.0Jb) and (c)]. 

If n > 2, the eigenstates are in general superpositions of 
Fock states. For n = 2 we get 



E2,i=K-6-2fx, |n = 2)i = |1)a|1) 



(22) 



'Uo-Ua 



2 V V 2 

\n = 2)± = c^^l |0)a|2)o + 4'2 |2)a|0)o 



-S] +P2_2^, 



J2) 



- 1 ± 




Uo-Ua- 26 



1/2 



It 



{Uo - Ua - 2Sy + 

Uq-Ua- 26 
{Un-UA-26f+4P\ 



1/2 
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In the case 9 



1/v^ and 



.(2) 

"2,± 



Tl/\/2. For 



/I > J7 — |P| the ground state is either |n = 2)_ or an arbitrary 
linear combination of \n = 2)i and |n = 2)„: If ga > 0, |n = 
2)_ is a ground state and the mean particle numbers = 
no = 1. If .g, < 0, E2,- = £2,1 = C/ - |P| - 2^ < £3^+, 
i.e., the ground state is doubly degenerate. Nevertheless, we 
have again n\ = tiq = 1. 

For n = 3 



■iUo + Ua 



K-25 



|n = 3),±=4'2|0)A|3)o + 4±|2)A|l 



K -S] + 3P2 _ 



,(3) 



„(3) _ 
— 



J3) _ 
^2,± — + 



1 ± 



3f/o -Ua-2K-2S 



It 



(3C/o - C/a - 2is: - 25)^ + 12P2 

3Uo -Ua-2K -26 
(3C/o - C/a - 2/^ - 2(5)^ + 12P2 



^'Sui = \- K - d 



±j(El^^+K-6] +3P2-3M 



|n = 3)„± = cf2 |1>a|2)o + 4± |3)a|0)o , 



„(3) 



1 ± 



Uo - SUa + 2K - 25 



1/2 




It 



(C/o - 3C/a + 2K- 25f + 12P2 

C/o - Wa + 2K-26 
(Uo - 3Ua + 2K - 2Sf + 12P2 



1/2 



In the case 6 — 0, the states are doubly degenerate and we 
have Esg^ = E^u- < E^g^ = £'3^+. If .ga < 0, the co- 
efficients are c^^l = cfj = 1/2, cfj = = ^3/2, and 
in spite of the degeneracy of the ground state we have ua = 

TiQ — 3/2. If ga > 0, the coefficients are Cq'^ — = \/3/2, 

c\_ = — 1/2, and the mean particle numbers ua, uq be- 
come undetermined. The degenerate ground states \n = 3)g_ 
and \n = 3)u_ occur in a domain for /i > 2{U — |P|), if 
5q < 0, and for ^ > 2U, if ga > 0. 



V. COLLECTIVE EXCITATIONS IN THE SUPERTLUID 
PHASE 



From the Hamiltonian (|9} we obtain the following Heisen- 
berg equations of motion: 



in,— a At 
at 



at 



-J A (aA,i-l + a-KA+l) + UAflAiO-M 



+ KfiQiCiAi - |P| a\^aQiao^ - /iflAi , (24) 

— — Jo (ao,i-l + flo^i+l) + J/o^-OiOoi 

+ KhAiCLQi - |P| SjiaAiftAi - ((5 + yu)aoi . 



These equations remain valid in higher-dimensional lattices as 
well with the only modification that the lattice indices aquire 
a vector character. 

In the mean-field approximation, the operators cia-i are re- 
placed by position independent c numbers a^' a^i ~ aa. One 
then finds a mean-field equation of the classical Hamiltonian 
form 



1/2 
(23) 



dt da*. 



(25) 



1/2 with the local classical mean-field Hamiltonian 



H = -2 (Jo|«oP + JAlaAp) + ^|«o|' + ^IttAl" 



9 |P| 

-XlaoCkAl^ Y [("a"o)^ + [oiIoiaY] 

-{S + fi) \aof - ^l\aA\'^ ■ 



(26) 



At least roughly, the correspond to the previously intro- 
duced order parameters ip^, but it should be noted that the 
present mean-field approximation is different from the one of 
the previous section rEa.(ll7H. where it was employed only to 
approximate the hopping term of the Hamiltonian, while now 
we use it equally in all terms of the Hamiltonian. The differ- 
ence of these approximations directly corresponds to the dif- 
ference of our goals in this and the preceeding section: there 
we aimed at a description of the Mott state, where the tp^r both 
vanish, and its eventual instabilty toward the formation of a 
superfluid state with an incipient (small) nonvanishing ipa-, 
while here our goal is to describe the excitations deep in the 
superfluid states, where the mean field is large and contains 
practically all particles. The price we have to pay for achiev- 
ing the latter goal is that the Mott state cannot be described by 
the modified form of the mean-field approximation. 

For 9 — Owe can use Eqs. il5i to simplify the local mean- 
field Hamiltonian, which can then be written as 

if = -(/i + 2J)n + ^n^ - ^ (^alao ± a*aAf , (27) 

where "+" and "— " correspond to 5^ < and ga > 0, respec- 
tively. 

In the equlibrium state, we set the time derivatives of the 
mean fields to zero. Then it is easy to verify that the phases 
of ao and aA must be equal, and may be set to zero for con- 
venience, i.e., we may, for the present purpose, simplify the 
mean-field ansatz 



(28) 
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with the normalization condition nA + ng = n. It is the small 
but nonzero value of |P| which aligns the phases of the a^ri 
for cr — 0, A. We are then left with the equations 

= (-2Ja + UaUa + Kuq - \P\no) ^/nX , (29) 
= (-2 Jo + C^o'^o + KnA - \P\nA - S) y/n^ . 

We consider the case ga > first. Then eas.(l29> in general 
have no positive solutions for no and ua with n = uq + ua- 
However, there is an exception for the special case = 
where, as we have discussed, the usual U(l) symmetry of par- 
ticle number conservation is enhanced to a U(l)xU(l) sym- 
metry. As it was also mentioned above, it is in any case rea- 
sonable to consider only small values of 9. Since the varia- 
tion of S with 9 is much faster than that of Uq, Ua, and P, 
for the semiquantitative analysis one can put Ua ^ Uq and 
K TnU^ + P. Then Eqs. (|29} take the form 

M\Am" = (-2 Ja + C/o") , (30) 
P^/n^ = (-2 Jo + ?7on - 5) . (31) 

Let us assume that ua ^ 0, uq ^ 0. Then we can di- 



vide Eqs. i3Q\ and (13 1> by -Jnp,y/riA, respectively, and we 
see immediately that Eqs. j30t,j31> are compartible only if 
J A = Jq + S/2, which does occur for = 0. Therefore in the 
degenerate case 9 = only the total density n = uq + ?ia is 
determined, 

^i = Un-2J , (32) 

and the actual values taken by uq and ua break spontaneously 
the SU(2) symmetry by fixing an arbitrary orientation of T in 
the (1, 3) plane. For 6* 7^ 0, either nA or no must vanish, i.e., 
T is then oriented parallel or antiparallel to the 3 direction, 
respectively. Since the mode has a lower energy, the solution 
of Eqs. j30t and i'ili is 

no = n , nA = , /i + (5 = Uqu — 2 Jo , (33) 

exactly as in the case of one-component BEC fT^flTil . 

In the case ga < 0, the situation is different and the solution 
of Eq. ( I29I 1 has the form 



nA = 



(C/o- JC+|P|)n-2(Jo- Ja)-<5 



Ua + Uo 
- (Jo + Ja) 



2{K-\P\) 



no 



(£/A-i^+|P|)n + 2(Jo- Ja) + ^ 
Ua + Uo -2{K- \P\) 



(34) 



2 (Jo - Ja) (t/o -Ua)-5 [Ua + \P\ - K) + 


UoUa-[K -\P\f 


n 


Ua + Uo-2{K-\P\) 





provided that 



{Uo -K+ \P\) n - 2 (Jo - Ja) - (5 > 



(35) 



otherwise the solution is given by Eq. 03> . For no and nA 
both nonvanishing the Bose gas has a nonvanishing hyperfine 
and magnetization density in the x direction proportional to 
^noUA- 

In the case ga < 0, 9 = 0, Eqs. ( I34> take a very simple 
form. 



no = nA = -, p={U -\P\)n-2J . (36) 



The excitations can be worked out employing the ansatz (3 



(37) 



Substituting Eq. ( I37t into Eq. MA\ we get the following lin- 



earized equations for the excitations: 

- Ja + ^KA+ij 

+ (2t/AnA + Kno - p) Aa» + {Uaua - \P\no) 
+ (if-2|P|)V^^Ao, + ifV^^Aj, , (38) 

- Ja ^Ao,i-i + Ao,i+i^ 

+ {2Uom + KuA - A* - (5) Aoi + (C/ono - |P|nA) AJ, 
+ {K- 2\P\) y/T^AA^ + K,/E^l\\^ . 

From Eqs. j38l l we obtain two coupled wave equations for the 

density fluctuations pai = y^n^ {^^cri + Aj^, 



ih—AAi 

at 



ih—Aoi 
at 



— Ja {P(yA-2 + PcrA+l) 



(39) 



1(0,1 



(0) 



The coefficients A)^^, ' are listed in the Appendix. The eigen- 
modes can be found by means of Fourier transformation, 

Pai = } , CakOk exp [I (kZi - UJkt)\ . (40) 



^ Cakbk exp [i {kzi - uJkt)] 

k 
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In the case 9 = 0, the two wave equations for the sum and the 
difference of the two densities are decoupled. The coefficients 
Cak are then independent of k. For ga < the two eigenmodes 
are given by 



huj[:^'> = ^£{k) [£{k) + 2{U -\P\)n] 



„(+) / J+) 



hujlr^ = £{k) 



2\P\n 



„(-) / J-) _ 



and for 5^ > by 

nJ^^ = ^/£{k) [£{k) + 2Un] 



„(+) /.(+) 



^£{k) [£{k) + 2Pn] , c^jC^/ci^ 



where 



£:(fc)=4Jsin^(|A 



(41) 



= 1 , (42) 



(43) 



is a dispersion relation of free particles in a lattice. The 
symmetric modes with the eigenfrequencies a;^^"* are quali- 
tatively the same in both cases {ga < and ga > 0). They 
are the Goldstone modes of the spontaneous breaking of the 
U(l) symmetry of particle number conservation, which corre- 
spond to condensate excitations. The antisymmetric mode is 
principally different: In the case ga > 0, it is an additional 
Goldstone mode [in which the components dT± transverse to 
the isospin vector in the (1, 3) plane oscillate] corresponding 
to the spontaneous breaking of the isospin symmetry, by the 
ground state; but for ga < the isospin is oriented along 
the 1 -direction by the interaction and the isospin symmetry 
around this axis is not spontaneously broken. In the anti- 
symmetric mode there is therefore a gap of size 2|P|n, and 
(cj[~^ — 2|P|ri) ^ k"^ for small k. This is again an isospin 
wave with ST2 and 6T3 propagating together. In the case 
ga < 0, the gap allows the energies of the symmetric and 
antisymmetric modes to cross at the points 



k 2 \P\^n 

— = ± arcsm — \ / ; — ;- 

kL n\2J{U-3\P\) 



(44) 



In the case 7^ 0, there are also two excitation modes. Let 
us consider the case ga < first. Typical k dependences for 
iiA ^ 0, no 7^ 0, are shown in Fig. ^| The two branches 
do not cross anymore and the coefficients c\k, cok depend 
strongly on k near the points of avoided crossings ( I44> . where 
the type of oscillations changes. If we increase 9 higher then 
a certain value the condition ( I35t is not fulfilled. According 
to Eq. ( I33> . Eqs. ( I38t become decoupled and the energies of 
the eigenmodes are given by 



Tiiuok = V£o{k) {£a{k) + 2Uon) , (45) 
huAk = [[£Aik)~iUo-K)n + 2{Jo-JA) + 6f 

-{Pnrf' , (46) 

where £o(A) is given by Eq. ( I43> with J replaced by Jo(A)- 
Since inequality i35\ does not hold and £a{Q) = 0, the A 
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FIG. 10: Spectrum of collective excitations hu>k and the ratio 



CAk/cok for **^Rb (ga < 0). q = CiI/AujrA 
n = 2. 



-4.5, 



0.7° 



mode has a gap at fc = 0. Moreover, [/q — (~ |-P|) is 
positive as well as all other summands in the square brackets 
in Eq. ( I46t . This means that for some nonvanishing values of k 
and large enough n the frequency uj\k can become imaginary, 
i.e., the system can be unstable. 

In the case ga > 0, 9 Q, the excitation modes are given by 
Eqs. i45i and ( I46> . and if — [/q (~ P) is positive. Therefore 
the frequencies cjAfc are always positive and there is a gap at 
k = 0. For 9 = this spectrum reduces to that given by 
Eq. ( I42I 1 as it should. 



VI. CONCLUSIONS 

We have investigated QPT of spin-1 bosons in an optical 
lattice created by the lin-6'-lin laser configuration. The laser 
beams create not only periodic potentials for atoms, but also 
couple the atomic ground states with magnetic quantum num- 
bers TO = ±L As a result there are only two lowest-energy 
Bloch bands (0 and A) which should be taken into account 
and the original three-component system is reduced to a two- 
component one. 

The interaction between the atoms is characterized by two 
scattering lengths, the symmetric gg and the much smaller but 
very important asymmetric ga- This setup could be realized 
with Rb or Na atoms, for which ga is negative (ferromag- 
netic interaction) and positive (antiferromagnetic interaction), 
respectively. It turned out that an appropriate change of the 
easily tunable parameters of the setup, the laser intensity, and 
the angle 9, permits a very rich scenario of phase transitions 
which we explored in detail in Section IV. The symmetries of 
the Hamiltonian are discussed and it is shown that the proper- 
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ties are essentially different for the atoms with ferromagnetic 
and antiferromagnetic interactions. Depending on the sign of 
the asymmetric scattering length, the isospin symmetry is or is 
not spontaneously broken in the superfluid phase. The corre- 
sponding collective modes, besides the always present Bogoli- 
ubov mode, are isospin waves and either gapless or gapped at 
k ^ 0. A surprising feature we found in this study is that the 
phase transition between Mott and superfluid phase can also 
be of first order in certain domains of parameter space, with 
associated domains where both phases can coexist, but, as we 
found, this is only possible for > 0. 



Expressions for Aqq', A^^^, aI^q , and A^^jl can be obtained 
from these equations making the replacement A — > 0, A, 



Acknowledgments 



This work has been supported by the SFB/TR 12 "Symme- 
tries and universalities in mesoscopic physics". 



APPENDIX: COEFFICIENTS IN EQ. 03 

The coefficients in Eq. ( I39> are given by 

a'^11 = -2JA{2UAnA + Kno-^) (A.l) 

- 2[Jo|P| - JA(i^- |F|)]nA (A.2) 

A^°l = 2Jl-4\P\{K-\P\)nonA 
+ [U^nA + iK+\P\)no-^i] 
x[3UAnA + iK-\P\)no-^i] (A.3) 

a'-^^ = 2nA {[U AHA + {K+\P\) no -^l]{K-\P\) 

-\P\ [SUono + {K- \P\) UA - (/i + <5)]} (A.4) 
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